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Abstract 



We investigate the influence of lunar-like satellites on the infrared orbital light 
curves of Earth-analog extra-solar planets. Such light curves will be obtained 
by NASA's Terrestrial Planet Finder (TPF) and ESA's Darwin missions as 
a consequence of repeat observations to confirm the companion status of a 
putative planet and to determine its orbit. We use an energy balance model 
to calculate disk-averaged infrared (bolometric) fluxes from planet-satellite 
systems over a full orbital period (one year). The satellites are assumed to 
lack an atmosphere, have a low thermal inertia like that of the Moon and 
span a range of plausible radii. The planets are assumed to have thermal 
and orbital properties that mimic those of the Earth while their obliquities 
and orbital longitudes of inferior conjunction remain free parameters. Even if 
the gross thermal properties of the planet can be independently constrained 
(e.g. via spectroscopy or visible-wavelength detection of specular glint from 
a surface ocean) only the largest (~ Mars-size) lunar-like satellites can be 
detected by light curve data from a TPF-like instrument (i.e. one that achieves 
a photometric signal-to- noise of 10-20 at infrared wavelengths). 

Non-detection of a lunar-like satellite can obfuscate the interpretation of a 
given system's infrared light curve so that it may resemble a single planet with 
high obliquity, different orbital longitude of vernal equinox relative to inferior 
conjunction and in some cases drastically different thermal characteristics. If 
the thermal properties of the planet are not independently established then 
the presence of a lunar-like satellite cannot be inferred from infrared data, thus 
demonstrating that photometric light curves alone can only be used for pre- 
liminary study and that the addition of spectroscopic data will be necessary 
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to properly characterize extra-solar Earth-like planets. 

Keywords: Planetary systems - planets and satellites: general - astrobiology 
- methods: data analysis 
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1. Introduction 



Planets with a minimum ma ss of 5-6 Ea r th masses have recently b een de 



tected around low mass stars (jUdry et al. 



2007 



Rivera et al 



2005|) and it 



seems li kely that obse r vator ies such as CoRoT or Kepler will detect yet smaller 



planets (iGillon et al. 



20051 ). Space-based observatories of the future will be 
capable of directly detecting Earth-sized planets around other stars. Pro- 
pos ed missions include a coronagraph operating at visible wavelengths (TPF- 



C) fjTraub et al 



20061 ) . and a 



infrared (TPF-I and Darwin) (jBeichman et al. 



arge-baseline in t erferometer operating in the 



2006 



Fridlund 



2000). One 



goal of such missions is to distinguish between planets that are Earth-like 
and can support life, and those that are decidedly less so (e.g., analogs to 
Mercury, Venus, or Mars). Several techniques have been proposed to carry 
out this classification. Spectroscopy can reveal the presence of atmospheric 
gases such as H 2 0, CH 4 an d O2, which are indicati ve of temperate conditions 
and/or biological activity (jpes Marais et all . I2OO2I ) . Photometry in reflected 
light can reveal diurnal (rotational) variability associated with ice, oceans, 
l and and vegetati on across the surface of a planet if no clouds are present 



( iFord et al. 



200ll ) . The specular "glint" from oceans might be detected as an 



increas e in the visible flux and polarization o 



angles ([Williams fc Gaidosl . 



2008 



McCuUough 



reflec te d light at lar ge phase 



2008|). 



Selsisl (12004 ) showed 



that orbital infrare d light curves could revea l general thermal properties of 



terrestrial planets. iGaidos fc Williams! (12004 . hereafter GW04) showed that 



diurnally averaged, orbital light curves at thermal infrared wavelengths con- 
tain information about the thermal properties of the planet's emitting layer 
(surface or clouds) and obliquity. Such light curves would be generated as a 



byproduct of repeated observations to confirm the companion status and or- 
bit of a putative planet and function as a first step towards characterization. 
One finding of GW04 was that oceans or a thick atmosphere damp seasonal 
variations in temperature, and that low or no variability in a planet's infrared 
light curve is indicative of the presence of oceans or a thick atmosphere. In 
conjunction with other characteristics, this is a signature of habitable surface 
conditions. These authors and the work presented here do not consider the 
effects of variable cloud cover. 

The disk-averaged infrared flux of an orbiting planet can vary as it presents 
different phases to a distant observer. This phenomenon has been observed 
for J upiter-mass extra - solar planets with s emi- major axes much less than 1 



AU ([Harrington et al. 



2006 



Cowan et al. 



2007 



Knutson et al 



2007J). The 



variation in flux from a planet depends on the diurnal pattern of outgoing 
infrared flux from the emitting surface (either the top of the atmosphere, 
if any is present, cloud layers, or the surface), which is controlled by the 
planet's thermal properties and day length. In general, significant day-night 
temperature differences will occur only if 



(1) 



OjJ < 



dl_\ 
df) 



where c is the heat capacity of the surface/atmosphere, u is the angular ro- 
tation rate, and the right hand side is the slope of the outgoing bolometric 
infrared flux at the emitting surface vs. temperature, T, evaluated at the mean 
surface temperature of the body, T. For the Earth: c = 8.34 x 10 7 J m~ 2 K _1 
and dl/dT = 1.58 W m~ 2 K _1 for T = 288 K (however the effective emitting 

temperature of the Earth is 255 K). Thus equation [1] does not hold for the 
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Earth: the Earth's day-night temperature variation is small and would be for 
any Earth-like planet with a rotation period much less than a year. This is 
primarily due to the high heat capacity of the ocean mixed layer which also 
moderates surface temperatures over landmasses and controls the outgoing 
infrared flux budget. 

The primordial rotation periods of terrestrial planets are thought to be a 
stochastic outcome of the final stages of formation by accretion o f plane tary 



2000). As 



embryos, and will be on the order of hours to days (ILissauer et all 
a consequence, the disk-averaged infrared flux from an Earth-like planet will 
only vary significantly along the orbit if the planet has a non-zero obliquity or 
eccentricity and hence seasons. This was explored in GW04. 

The Earth's Moon, lacking an atmosphere or oceans and having a lunar day 
29.5 times longer than the Earth, experiences a much larger diurnal surface 
temperature variation. Absence of recent geologic activity on the Moon has 
allowed a regolith of impact ejecta to accum ulate. This ma t erial is optically 



dark [the average lunar Bond albedo is 0.07 (ILane fc Irvind . 119731 ) compared 



to the Earth's 0.31] and has a relatively low heat capacit y [that of the Moon 
is 4 x 10 4 J m" 2 K" 1 at 29.5 days or 0.1% of the Earth, (jMuller fc Lagerros . 



19981 1]. As a result, the inequality of Equation[T]is satisfied and thus the Moon 
makes a significant or even dominant contribution (depending upon viewing 
geometry) to the variable component of the infrared flux from the Earth- Moon 
system. 

The Moon is thought to have accreted from a ci rcumterres t rial d isk of 



1986|). The 



ejecta generated by the impact of a Mars-sized body (lHartmannl . 
high ambient temperatures and low gravity in the transient disk explain the 
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Moon's lack of volatiles (IPritchard Stevenson 



20001 ). Current scenarios 



for the final stages of t errestrial planet formation include such giant impacts 



(jCanup fc Agnor . 



19981) and the results o 



they are not rare (jlda. Canup fc Stewart 



nume rical simulations suggest that 



19971 ). Thus large satellites lack- 



ing atmospheres or oceans may be common around extra-solar rocky planets. 
Like the Moon, these satellites would have originally formed closer to their 
parent planets and their rotations would haye quic kly synchronized to their 



orbits ( Gladman et al. 



1996; 



Canup fc Agnor 



19981 ) . As a result their diurnal 



temperature variation could be significant. We note that Moon-size satellites 
could retain an atmosphere against gravitational escape over Gyr time-scales 
if one was originally present. 

Satellites around extra-solar planets will be unresolved by even the most 
ambitious planet-finding mission: the angular separation of the Earth and 
Moon at a distance of 10 pc is 0.25 mas. However, a large satellite might 
reveal itself by a significant variation in the total (bolometric) flux from the 
system. Such an interpretation requires independent knowledge of the ther- 



mal and rotational properties of the parent 



using spectroscopy ( 



Pes Marais et al 



data (IWilliams fc Gaidos 



2008 



2002 



planet, whi ch can be established 



Sclsis. 



McCullough 



20041 ) or optical light curve 



20081 ). If establishment of these 



gross thermal properties leads to the expectation that infrared flux variation 
would be small (Equation [1]) then observation of significant variation could be 
attributed to the presence of a large satellite. In the absence of such auxiliary 
information, however, the satellite contribution may result in an assignment 
of erroneous thermal properties to the planet. 
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We present calculations of infrared light curves of an Earth-like planet with 
a Moon-like satellite. The terms Earth-like or Earth-analog refer to specific 
thermal and orbital properties that represent those of the Earth (see The 
albedo, heat capacity and orbital period of the satellite are set to that of the 
Moon (see above) while its radius is allowed to vary. 

In §2] we describe the details of the analytical energy balance model (EBM) 
used in these calculations. We give an illustrative calculation in §31 We then 
determine the minimum radius of lunar-like satellite that can be detected at 
infrared wavelengths around an Earth-analog planet O- In §5]we describe the 
biases in planetary orbital properties that can be introduced by an undetected 
lunar-like satellite. In §6] we describe the effects that low and high altitude 
clouds would have on our calculations and in JT] we discuss the implications of 
our results. 



2. Model 

Our calculations are based on the infrared orbital light curve model of 
GW04. These authors employ a linearized, analytic EBM to calculate the 
infrared flux emitted by a planet. This model assumes a single, uniform plan- 
etary albedo and parameterizes the thermal inertia and meridional heat trans- 
port across a planet's surface. The thermal effect of clo uds is accounted for 



by su btracting a correction term from the outgoing flux (ICaldeira fc Kasting , 



19921 ). The time-dependent surface temperature distribution is described by 
a combination of Legendre polynomials and a Fourier series that are solutions 
to a diffusion equation with periodic temporal boundary conditions in a spher- 
ical coordinate system. The disk- and diurnally-averaged infrared flux for a 
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prescribed viewing geometry is calculated along an entire orbit. As long as 
Eqn. [T]is satisfied, then diurnally averaging the infrared flux justifies the use 
of a single, average planetary albedo. 

We consider only one set of Earth-like planetary parameters. Although such 
properties will undoubtedly vary amongst extra-solar planets, those planets 
with thermal properties similar to the Earth will be the most compelling tar- 
gets of investigation. The thermal inertia of the surface (8.34 x 10 7 J m~ 2 K _1 ) 
and heat diffusion coefficient (0.38 W m -2 K _1 ) were chosen so that with an 
albedo (A = 0.3055), orbital semi-major axis (a = 1 AU), and eccentricity 
(e = 0.0167) of Earth, the model reproduces the meridional surface tempera- 
ture distribution of the Earth as well as the seasonal temperature variation at 
several latitudes (GW04). 

The orbital properties (a, e and i) of a real planet can be determined by 
imaging or astrometry, but the thermal properties of the planet may not be 
uniquely determined by independent means. We discuss this scenario in §0 
Under the conditions of known orbital and thermal properties, the light curve 
is a function of the planet's obliquity (<5 ), the orbital longitude of inferior 
conjunction relative to the spring equinox (L ), and the orbital longitude of the 
apastron (L ap ). If the orbit of the planet is nearly circular then the longitude 
of apoastron L ap (fixed here to 180°) is unimportant. 

Ocean and atmospheric circulation and the thickness of an ocean's mixed 
layer may differ for Earth-like planets with obliquities that are significantly 
larger than 23.5°. Thus, the actual light curves of such planets would differ 
from those calculated with Earth-like thermal properties. We examined this 
effect by comparing general circulation model (GCM) runs for 5 = 85° to 
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our EBM calculatio ns (Figure ED). The GCM used was the three-dimensional 



GENESIS 2 model dWilliams fc Pollard 



20031 ) and the calculations were per- 



formed with L = 120°. The amplitude and general shape of the GCM light 
curves are nearly identical to those of the EBM. However, we find that the 
phase of the GCM light curves tend to lead those of the EBM by ~ 35°, with 
the greatest differences occurring when the geometry of the system is such that 
the poles of the planet are pointed towards the observer (e. g. high obliquity 
and high inclination). This is likely due to the inclusion of polar sea ice in the 
GCM. 

The exact origin of the offset in phase is uncertain, but we suspect that it 
involves effects of seasonal changes in cloud cover, sea ice or ocean circulation 
not included in the EBM. In §Hwe show that artificially adjusting the phase of 
the EBM light curve (for i = 60°) to better match the GCM calculations only 
slightly increases the probability of satellite detection. In £0, we find that this 
phase lag does not affect the conclusion that an Earth-like planet with a lunar- 
like satellite produces light curves resembling those of an isolated planet with 
high obliquity (i.e. one with large amplitude, see GW04 for for a comparison 
of light curves from Earth-like planets with high and low obliquities). We 
are interested only in estimating the detectability of satellites and their gross 
effect on the interpretation of infrared light curves, rather than on detailed 
inferences about the climates of planets themselves and thus use the EBM to 
efficiently calculate light curves over a range of obliquity values. 

Our formalism for calculating the outgoing infrared flux from a lunar-like 
satellite is also based upon the analysis in GW04. In the absence of an atmo- 
sphere or oceans, the energy-balance equation governing the temperature (T) 

10 



at a given point on the surface of a satellite with no latitudinal heat transfer 
is 



(2) 



c- 



dT 
dt 



S-(1-A)-I(T), 



where time is denoted by t, incident stellar flux by S, albedo by A, and out- 
going infrared flux I(T). S and I(T) are calculated as functions of longitude 
and latitude on the surface of the satellite taking into account projection ef- 
fects. To analytically solve this equation three assumptions are made. First 
c is assumed constant in time and across the satellite surface. Second, the 
temperature dependence of the outgoing infrared flux is approximated as a 
linearized blackbody: 



(3) 



I(T) = I(T) ■ (l + 4(T - T) /T) . 



This follows the approach of classical energy balance models ( jNorth et al. 



19811 ). Finally, we assume tidally-locke d, synchronous rotation as is expectec 



for la rge, collisionally-formed satellites (IGladman et al. 



1996 



Canup fc Agnor 



19981 ). These assumptions allow Fourier series solutions to Equation EJ 

N 

(4) T(6, 1) = T {6) + J>n(0) cos(n£) + b n (6) sm(n£)], 

n=l 

where T is the mean temperature for a given latitude 9 and i is the longi- 
tude on the surface of the satellite. N is set to 10 as numerical tests show 
that larger values do not significantly change the final light curve. Substitu- 
tion of Equations [3] and H] into Equation [2] yields expressions for the Fourier 



n 



coefficients: 



(5) mid = <{ ^- u :° 

cr + (B + a) 

(6) "»(») = ^T^V 

cr + [B + a) 2 

where a' n , b' n , B and a are 

(7J a n = cos{ni) at 

Jo ^ 



o * 



4/(T) 

(9) B = -U. 



(10) a = ^v^. 

The infrared emission from the surface is calculated using Equation [3] and the 
total signal is determined by geometric projection of the hemisphere facing the 
observer. 

For simplicity we assume that the orbit of the satellite is coplanar with that 
of the planet's orbit around the star and, because of synchronous rotation, 
the satellite has zero obliquity. Thus any variation in outgoing flux from 
the satellite is due to its finite heat capacity. The disk-averaged flux of the 
satellite is then independent of the location on its orbit around the planet 
and depends only on the geometric angle described by the star, satellite, and 
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distant observer (Figure [2]). Variability in the satellite signal is due to its 
observed phase, which changes with the orbital period of the planet. In §H 
we consider satellites that differ in size (but not surface properties) from the 
Moon. 

A satellite larger than the Moon will retain heat for a longer time and be 
more likely to have active volcanism. This could make the satellite darker, as 
in the case of the lunar Mare. However, with an average albedo of 0.07 the 
Moon is already quite dark. Fresh basalt from active volcanism on the surface 
of a larger satellite would have little effect on its light curve. 

3. Example Light Curves and Observations 

Figure [2] illustrates how an Earth-Moon analog would appear at five evenly- 
spaced points in the system's orbit. Figure [3] plots the infrared light curves 
produced by this system. The bottom panel displays the disk-averaged flux, 
while the top shows light curves normalized to their respective means. This 
normalization (which is employed for all subsequent analysis) removes the 
radius of the planet as a degree of freedom in the model. The calculations were 
performed assuming an Earth "twin" (5 = 23.45°) with a satellite of radius, 
orbital period and albedo equal to that of the Moon (0.273 R e , 29.5 days and 
0.073 respectively). For these and all simulations the coplanar orbits of the 
planet and satellite are inclined by 60° with respect to the plane of the sky 
(the median value of an isotropic distribution). The dotted line in both panels 
is the contribution from the planet alone. The peak-to-peak amplitude of the 
satellite's signal (dash-dots in the bottom panel) is 55 W m~ 2 , whereas the 
planet's flux alone varies by only 4 W m -2 . Because the thermal inertia of 
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the satellite is low, it displays a larger relative infrared flux variation than 
that of the planet (Equation [TJ). However, because the satellite is smaller, the 
majority of the average flux originates from the planet. 

The assignment of erroneous properties to a planet with an undetected satel- 
lite is illustrated in the top panel of Figure [3j where the dashed line is the 
best-fit planet-only model to the planet+satellite observations. This model is 
of a planet with Earth-like thermal properties and 5q = 75° and L = 90°, i.e. 
quite different than the input values. The phase and amplitude of the light 
curve produced by these models are not independent and each depends on the 
obliquity and orbital longitude of inferior conjunction. Adjusting the orbital 
longitude of inferior conjunction could produce better agreement with phase, 
but would unacceptably decrease the amplitude of the light curve. In general 
large light curve amplitudes like those produced by Earth-like planets with 
Moon-like satellites can only be mimicked by single planets with very high 
obliquities and orbital geometries where the northern or southern hemisphere 
is pointed towards the observer. 

An example of an observation scheme of five evenly-spaced measurements is 
shown in the top panel of Figure [3J In practice a minimum of three observa- 
tions are required to co nfirm planetary status and reject background sources 



( IBeichman et al 



20061 ). The TPF mission would conduct a minimum of 3-5 



observations on each star during the first two years of a five-year mission. The 



rem aining time wou 



ets ( IBeichman et al. 



d be sp ent on spectroscopic follow-up of a few dozen plan- 



20061 ). Spectroscopic observations could, in principle, be 
split into multiple integrations, however these would preferentially occur when 
the planet was near maximum elongation from the star, thus maximizing the 
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S/N. Both a nulling interferometer (TPF-I, Darwin) or a coronographic im- 
ager ( TPF-C ) will obscure planets along some parts of orbits in the habitable 



zone (IBrown 



20041 ). We thus consider two observing scenarios, one consisting 



of five observations at equal longitudinal intervals around the orbit, and a sec- 
ond consisting of 14 observation points restricted to half of the orbit furthest 
from the star. For example, 14 x 2-day integrations of 36 high-priority targets 
might be obtained in three years . Targ ets in the habitable zone of solar-mass 



stars (0.9-1.3 AU, 



Kasting et al. 



(119931 )) will have orbital periods of 300-550 



days and thus the accessible part of an orbit can be completely observed. For 
the 5-point "confirmation" observation scheme we assume a S/N ratio of 10 
per observation, the median value of the S/N amongst all 234 stars for which 
S/N > 5 is achievable for an Earth-sized planet in a 24-hour integration time. 
For the 14-point "characterization" scheme we assume a S/N of 20, which will 
be the case for Earth-s ized planets around the nearest 20% of the target stars 



( jBeichman et al. 



20061 ) . 



4. Satellite Detection Limits 

We determine the minimum size of a Moon-like satellite that can be detected 
around an Earth-like planet whose gross thermal properties are outlined in §5J 
We presume that the inclination of the planet's orbit with respect to the plane 
of the sky is independently measured by astrometry and use a single value of 
60°. For these calculations the orbital and thermal parameters of the planet 
and the orbital period, albedo and thermal inertia of the satellite (assumed 
to be equal to that of the Moon) are held fixed, while 5q, Lq and the satellite 
radius (R s ) are allowed to vary. With the exception of cases with extremely 
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high planetary obliquity where most of the planetary signal is at twice the 
orbital period (GW04), the signal from the planet and satellite will have the 
same period. There are then three unknowns (So, L and R s ) but only two 
measurable quantities; the amplitude and phase of the orbital signal. Thus 
it is not possible to uniquely disentangle the planetary and satellite signals. 
Instead, we define a satellite "detection" as the case where the observations 
cannot be accounted for by a planet-only model. 

Our detection analysis is as follows: We generate an array of planetary light 
curves over the full ranges of So [0-90°] and Lq [0-360°]. The satellite light 
curve for a given R s is calculated and added to each planet light curve in the 
array. Each total light curve is sampled at N specified points according to 
either of the observing schemes described in Random noise with a given 
RMS is added to these measurements. We then perform an exhaustive search 
of planet-only light curves to find the minimum x 2_ fit to the measurements. 
The analysis is repeated for different values of R s . We interpret the confidence 
level C associated with the value of x 2 an d the number of degrees of freedom 
as the probability that the deviation from the planet-only model is due to 
the presence of a satellite. This is because 1 — C is the probability that 
measurements of the light curve of the planet alone would result in a fit with 
a x 2 larger than the observed value, i.e. a false positive. We set the effective 
number of degrees of freedom to the minimum x 2 value associated with a best 
fit to N observations of the planet-only light curve. Although this minimum x 2 
value is sensitive to the random noise characteristics of a given sampling, we 
find that a more robust calculation of the true x 2 minimum (i.e. averaging over 
a large number of random noise patterns) produces results that are in close 
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agreement with this single noise characteristic approximation. Our minimum 
X 2 approximation produces an increase of nearly an order of magnitude in 
computational efficiency over the robust method. 

We average C over L Q , which cannot be independently determined for any 
system, but will have a uniform probability distribution. In addition, C is 
averaged over all possible phases of the measurement scenarios, i.e., the longi- 
tudes at which the planet is observed are shifted by 2 degree increments over 
the orbital phase range 2ir/N. This produces an average probability of detec- 
tion (C) for a given R s and Sq. The planet's obliquity cannot be independently 



determined from the light curve o 



to have an isotropic distribution (lAtobe fc Ida 



the planet ^- s atellit e, nor can it be assumed 



20071 ) 



In Figures H] and [5] we plot (C) vs. obliquity for several values of satellite 
radius (in units of planetary radius). For an Earth-sized planet, this range 
of radii corresponds to Vesta- to Mars-sized bodies. The two figures corre- 
spond to the "confirmation" (5 evenly-spaced observations at S/N = 10) and 
"characterization" (14 observations at S/N = 20 distributed around 50% of 
the planet's orbit) observation scenarios described in £j3l As the radius of 
the satellite increases the observations become increasingly inconsistent with 
a planet-only light curve, thus the probability of satellite detection increases. 
For satellites smaller than 0.33 planetary radii the probability of detection 
is lower in the "characterization" observing scheme. This is an effect of the 
observations being distributed around only 50% of the orbit. The incomplete 
phase sampling of this scheme does not capture the peak of the satellite flux 
which occurs at the longitude of superior conjunction. However for large satel- 
lites (R s > 0.33 planetary radii) the amplitude of the net light curve becomes 

17 



so great that the satellite is detected even without complete phase coverage. 
The scatter in these probability curves is due to the stochastic noise added 
to each of the sample measurements. For both observation schemes < C > 
is very weakly dependent on obliquity. These results show that a Moon-like 
satellite (0.27 Earth radii) would only be detectable with ~30% confidence 
by either observation scheme. For the "confirmation" and "characterization" 
observing schemes a satellite would have to be 0.5 and 0.38 planetary radii 
respectively to be detected with 90% confidence. This corresponds to approx- 
imately Mars-sized satellites in orbit around an Earth. 

In Figure [6] we consider how these results would change if a full three- 
dimensional climate model were used to generate the planetary light curves. 
As previously stated, our EBM produces light curves that lag by ~ 35° rel- 
ative to those of the GENESIS 2 GCM for obliquities of 23.5° and 85°. To 
mimic the results of the GCM we offset the phase of the EBM light curves 
by 35° and repeat the analysis of Figure [51 It would be computationally pro- 
hibitive to generate GCM light curves for the full range of obliquities that are 
included in this analysis, thus we approximate the GCM by applying this off- 
set. We find that the probability of detection actually increases slightly with 
the phase-adjusted pseudo-GCM light curves (Figure E]). Thus, we conclude 
that our EBM results are conservative estimates for the probability of satellite 
detection. 

5. Errors Introduced by an Undetected Satellites 

As we showed in Figure El the application of planet-only models to a set 
of observations can result in mischaracterization of the Earth-like planet if a 
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large satellite is present. Even if the gross thermal properties of the planet 
have been independently established, this will still produce erroneous values 
of Sq and Lq. We describe this effect by recording the "true" (<5 , L ) pair 
for a given planet+satellite light curve and the (5 , L ) pair of the best fit 
planet-only light curve as determined by a x 2 analysis. In some instances, a 
planet+satellite light curve can be fit with equally low \ 2 by more than one 
planet-only light curve. For these cases, we choose the (8 , L ) pairs that 
are closest to the true value. Figure [7J plots the direction and proportional 
magnitude of the error introduced by satellite confusion (The length of the 
vectors have been reduced for clarity). These simulations were run for an 
Earth-like planet with a satellite of radius equal to that of the Moon and (for 
clarity) no intrinsic noise added to the sample points. In nearly all cases, 
the presence of a lunar-like satellite makes an Earth-analog planet appear to 
have high (> 80°) obliquity. Shifting the planetary light curve phases by 35° 
to make them resemble the GCM results does not alter the tendency of the 
best-fit solutions towards high 5q. 

The vectors in Figure [7] converge towards two values of L$ (90° and 270°) 
because the planet-only light curve is in phase with that of the satellite at 
these L . For high obliquity planets, L Q values of 90° or 270° translate to a 
planetary configuration in which the southern or northern hemispheres (re- 
spectively) face the observer. At high Sq, these are the only geometries that 
produce planet-only light curves with large amplitudes that can reasonably 
fit a planet+satellite system. Adjusting for the phase difference between the 
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EBM and the GCM does effect these general trends in L . Instead of converg- 
ing to Lq values of 90° and 270°, the GCM vector field reaches convergence at 
Lq values that are shifted by 35°. 



6. The Effect of Clouds 



Clouds play an important role in the energy balance and climate of the 
Earth by reflecting sunlight and scattering and trapping long-wavelength radi- 
ation; they would presumably do so for Earth-like extra-solar planets as well. 
Althou gh clouds represe nt a mean 20 W m~ 2 (8%) gain in radiation for the 



Earth (jHartmannl . 



19941 ) we are concerned here only with their seasonal affect 
on the energy budget and disk-integrated outgoing radiation. These seasonal 
effects will be most prominent when one hemisphere is presented to the ob- 
server. On the Earth, low clouds (stratocumulus) produ ce a decrease in n et 



radiation, while high clouds (cirrus) produce an increase (jHartmannl . 1 19941 ). 

We estimate the thermal effect of seasonal variability in stratocumulus 
clouds on the light curve of an Earth twin. We assume that these low-altitude 
clouds radiate at the same temperature as a clear atmosphere and surface, i.e. 
they produ c e no a dditional greenhouse eff ect. We use a cloud albedo of 0.7 



(Hartmann 



199J). 



Comiso fc Stock! (120011 ) estimated variation in cloud cover 



over the open ocean around Antarctica to vary by only ±1% around a mean 
of 91%. Seasonal variation of ±10% around a mean of 80% was observed over 



the North Atlantic (IMassons et al. 



19981 ). We adopt a 20% seasonal variation 



in the mean as a reasonable bracketing value. To produce the global mean 
Earth albedo of 0.31, a dark ocean (A = 0.07) must be covered with 37% 
clouds. Although the average cloud cover on Earth is ~ 60%, only half of 
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those are low-altitude stratocumulus (IMinnis et al 



2002), a value in reason- 



able agreement with the calculated fraction of 37%. A ± 20% fluctuation in 
the mean coverage produces an albedo variation of 0.047. We assume that 
the albedo variation is uniformly distributed over each (northern/southern) 
hemisphere and that it varies sinusoidally in phase with the summer solstice. 
(This obviously produces a non-physical discontinuity at the equator which 
is unimportant for the purposes of estimating the magnitude of the effect of 
clouds). We examined the light curves of two cases: one in which an Earth 
twin is observed at a moderate inclination (i = 60°) and the other in which 
a high-obliquity Earth is observed on an edge-on orbit (i = 90°) so that one 
hemisphere is seen nearly pole-on (results not shown). In both cases the annual 
mean of the disk-averaged flux is slightly lower in the presence of low altitude 
clouds, but the amplitude and phase of the variation is essentially unchanged. 
Of course, pathological deviations from terrestrial patterns of cloudiness are 
possible on planets not quite like the Earth, but we have already shown with a 
GCM model (that includes parameterized cloudiness) that in at least the high 
obliquity regime, our conclusions are not significantly impacted. 

For an Earth-like planet, high altitude clouds will cause a greenhouse effect 
whose net effect is to offset any decrease in infrared emission caused by their 
high albedo. Variation in the fraction of high-altitude cloud cover on time 
scales of less than ~ 1 day will be averaged over during typical integrations. 
We argue that variation in high cloud cover over longer time scales will be 
insignificant compared to estimated noise characteristics for a typical observa- 
tion. High altitude clouds with a temperature of 210 K would emit ~ 40% less 
radiation than the surface. With typical high cloud cover fractions of ~ 30% 
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(IMinnis et al.l . 120021 ) and again assuming ± 20% fluctuations in coverage, we 
calculate that fluctuations in outgoing infrared flux due to high cloud variabil- 
ity will be on the order of a few percent (40% of surface flux x 30% coverage 
x 20% variability), which would be unresolved by observations with an opti- 
mistic S/N of 20. In a more extreme case (e.g. larger variability, greater mean 
surface coverage or lower emitting temperatures) the variability in a light curve 
due to high altitude clouds will act to further confuse its interpretation. 



7. Discussion 

The simulations presented here show that time-series infrared photometry 
by a TPF- or Darwin-like observatory would reveal only the very largest lunar- 
like satellites (Mars-sized) around Earth-analog planets, and then only if these 
Earth-like properties, i.e., the presence of oceans and/or a substantial atmo- 
sphere, have been established by independent means, e.g. spectroscopy or op- 
tical photometry. This conclusion holds for a wide range of planetary obliquity 
(Figs. EEf5f6j) . assuming that the approximations of the EBM do not grossly 
misrepresent the infrared light curve of a high-obliquity Earth-like planet (Fig. 

nj. 

When interpreting infrared light curves, the presence of an undetected lunar- 
like satellite can suggest erroneous values of planetary obliquity and longitude 
of inferior conjunction. In the case of a planet with high thermal inertia, 
inferred values of 5q near 90° and L values within 35° of 90° or 270° may 
indicate the presence of a lunar-like satellite. This result is based on the as- 
sumption that the satellite and planetary orbits are c oplan ar, which may be 



not be the case for high obliquity planets. 
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Kinoshital (119931 ) showed that the 



orbit of a satellite will stay in the equatorial plane of its host planet if the 
secular rate of change of the planet's obliquity is slower than the precessional 
speed of the satellite orbital plane. Thus satellites around planets that ex- 
perience d rapid changes in obli quity [possibly by collisions as in the case of 



Uranus (IParisi fc Brunini 



19971 )] would stay in their coplanar orbits. 



If a satellite's orbit is non-coplanar then its rotation axis will be tilted with 
respect to the plane of the planet's orbit. This effectively causes a non-zero 
obliquity for the satellite which will modify the amplitude of the satellite's 
light curve but will not change its period. If this non-coplanar orbit precesses 
then the signal will change over a time-scale of many orbital periods. 

If the ther mal properties of a plane t are not independently established via 



spectroscopy (IDes Marais et al. 



20021 ). visible- wavelength detection of glint 



from an ocean or significa nt polarization of visible reflectance (IWilliams fc Gaidos . 



2008 



McCullough 



20081 ). then the flux from an unresolved lunar-like satellite 



can induce serious errors. If the measurements are modeled with the several 
free parameters (e.g. 5q, Lq, A, c and efficiency of meridional heat trans- 
port) then a set of planet+satellite measurements can be satisfactorily fit by a 
planet-only light curve. For instance, the peak-to-peak light curve amplitude 
from a system with a large, unresolved lunar-like satellite around an Earth- 
analog planet can be fit by a planet with low thermal inertia and drastically 
different So and Lq, implying a planet more akin to Mars than Earth. Such an 
erroneous inference would impact the determination of the frequency of hab- 
itable planets. This reinforces the need for multiple wavelength observations 
including spectroscopy and phot ometry to di s ambiguate the characterizat ion 



of extra-solar terrestrial planets (ITraub et al. 



2006 



Beichman et al. 



20061 ) 



2:5 



Where, optimistically, the thermal properties of a planet are known and its 
satellite is Mars-sized, the existence of a satellite may be inferred from in- 
frared data. Such a discovery would provide information about the collisional 
and kinematic evolution of the parent planet. In addition, a large satellite 
could be a potential indicator of habitability, as the pres ence of the Moon i s 



known to stabilize the obliquity and climate of the Earth (ILaskar et al. 



1993|). 
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Figure 1 . Calculated EBM (solid) and GCM (diamonds) light 
curves of a high-obliquity (So = 85°) planet with Earth-like ther- 
mal properties. Light curves are shown for three values of incli- 
nation. L = 120° (defined with respect to the vernal equinox) 
is used for both models. The two models are in close agreement 
regarding the amplitude of the planetary signal, however the 
phase of the EBM calculations tend to lag by ~ 35°. 
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Figure 2. Planet and satellite at five uniformly-spaced orbital 
longitudes. The vernal equinox corresponds to L = and in- 
ferior conjunction to L — L . The disk-averaged flux from the 
Earth-like planet, which has a high thermal inertia, varies only 
with seasonal surface temperature differences between the two 
hemispheres. The disk-average flux from the Moon-like satellite, 
which has a low thermal inertia, depends on its observed phase. 
The observed phase of the satellite depends only on its orbit 
around the star, not around the planet. 
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Figure 3. Simulated light curves and observations of an Earth- 
Moon "twin" . The bottom panel shows disk-averaged flux while 
the top displays the normalized light curves. In both panels the 
dotted curve is the planetary seasonal flux and the solid curve 
is the "true" planet+satellite light curve. The dash-dot curve in 
the bottom panel is the satellite contribution. In the top panel 
the sample points are indicated by X's with one-sigma error 
bars for a S/N = 10. The best fit, planet-only light curve to the 
sample points is represented by the dashed line. The top panel 
light curves have been normalized by their mean to remove any 
dependence on the radius of the planet. 
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Figure 4. The L averaged probability of satellite detection 
as a function of planetary obliquity for a range of satellite radii 
(0.16, 0.22, 0.27, 0.33, 0.38 and 0.44 planetary radii). The radius 
of the Moon is 0.27 R e . These simulations were performed with 
5 evenly-spaced observations and a S/N of 10. As the radius of 
the satellite increases it becomes increasingly difficult to explain 
the sample measurements with a planet-only light curve, thus 
the probability of detection increases. 
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Figure 5. Similar to Figure HI except for the "characteriza- 
tion" observing scheme: 14 sample points with a S/N of 20, 
distributed around 50% of the planet's orbit. Although the S/N 
is higher than the observations of Figure HI the distribution of 
sample points around only 50% of the orbit makes it difficult to 
detect satellites of small radii (< 0.33 planetary radii). 
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Figure 6. Similar to Figure [5l except that a 35° phase shift 
has been applied to all EBM planetary light curves so that they 
agree with the GENESIS 2 GCM. Using these pseudo-GCM light 
curves the detection probability actually increases relative to the 
EBM case (Figure |5]). This suggests that the EBM places a lower 
limit on the probability of satellite detection. 
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Figure 7. Error induced in the light curve analysis by an un- 
detected satellite. The vectors point from the true (So, L ) pair 
that is used to generate the planet+satellite light curve, towards 
the best fit planet-only (S , L ) pair. For clarity, the length of 
the vectors have been scaled to the amount of induced error. 
Their endpoints do not actually fall on the best-fit (So, L ) val- 
ues. Nearly all best fit values cluster at high obliquity and L 
values of ~ 90° or ~ 270°. 
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